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Abstract. We present a detailed analysis of the dynamical response of ultra- 
cold bosonic atoms in a one-dimensional optical lattice subjected to a periodic 
modulation of the lattice depth. Following the experimental realization by Stoferle 
et al [Phys. Rev. Lett. 92, 130403 (2004)] we study the excitation spectrum 
of the system as revealed by the response of the total energy as a function of 
' ^ ' the modulation frequency f2. By using the Time Evolving Block Decimation 

C ' algorithm, we are able to simulate one-dimensional systems comparable in size to 

O ' those in the experiment, with harmonic trapping and across many lattice depths 

O , ranging from the Mott-insulator to the superfluid regime. Our results produce 

many of the features seen in the experiment, namely a broad response in the 
superfluid regime, and narrow discrete resonances in the Mott-insulator regime. 
^^ ' We identify several signatures of the superfiuid-Mott insulator transition that are 

*^ ' manifested in the spectrum as it evolves from one limit to the other. 
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1. Introduction 

The realization of ultracold atoms confined in optical lattices has made a large range 
of fundamental equilibrium and dynamical phenomena of degenerate quantum gases 
experimentally accessible. The success of this approach stems from the fact that, in 
contrast to analogous condensed matter systems, optical lattices form a defect free 
lattice potential which can trap a dense cloud of atoms with long decoherence times 
and can be controlled rapidly with a great deal of flexibility. This has already enabled 
a seminal demonstration of the superfluid (SF) to Mott insulator (MI) transition 
by Greiner et al jjj that was predicted to occur for a clean realization of the Bose- 
Hubbard model (BHM) |21 El ■ More recently, other themes of ultracold-atom research 
have been explored experimentally such as the purely one-dimensional (ID) Tonks- 
Girardeau limit 01 El , the characterization of the SF-MI transition via the excitation 
spectrum through ID - 3D dimensionality crossover jHllZlj s^nd impurity effects caused 
by Bose- Fermi mixtures El ■ 

Here we focus on features of the excitation spectrum from reference which were 
revealed for a ID BHM by periodic modulation of the lattice depth. The experimental 
accessibility of the excitation spectrum provides a rich source of additional information 
that can be compared with well studied quantities such as the dynamic structure factor 
|l()l ITT] . In addition to this the experiment demonstrates the transition through the 
evolution of the spectrum from discrete sharp resonances in the MI regime to a broad 
continuum of excitations in the SF regime. Changes in the structure of the excitation 
spectrum provide important evidence for the transition beyond the loss and revival 
of phase coherence when ramping the lattice 0, and can also be used to diagnose 
the temperature of the system ^1 El ■ The ID system is of particular interest for 
several reasons. Firstly, quantum fluctuations are expected to play a strong role 
there jjll, and this is indeed found to be the case with the critical ratio of the on- 
site interaction to kinetic energy {U/ J)c, identified by the appearance of the discrete 
structure in the spectrum, being lower than that predicted by mean field theory |^. 
Secondly, the behavior found in the ID experiment for the SF regime, specifically a 
large and broad non-zero response, most strikingly departs from standard theoretical 
predictions. Specifically, linear response using Bogoliubov theory for a Bose-Einstein 
condensate (BEG) in a shallow ID optical lattice predicts that lattice modulation 
cannot excite the gas in the SF regime due to the phonon nature of the excitation 
spectrum 15 . Since the quantum depletion of the SF in the experiment was significant 
(~ 50%), it has been suggested that this was responsible for the response jl£}. Indeed, 
it has been shown that only a small amount of seed depletion is required for non- 
linear effects like the parametric amplification of Bogoliubov modes to reproduce 
the SF response |17[ I18| . More recently still, the use of the sine-Gordon model 
and bosonization method has demonstrated that linear response is non-zero at low 
frequencies ^Hl- Lastly, the study of the ID system permits the use of quasi-exact 
numerical methods, such as the Time Evolving Block Decimation (TEBD) algorithm, 
where the fully time-dependent dynamical evolution of the system can be computed 
efficiently for systems of equivalent size to those in the experiment. 

In addition to the many-body physics perspective, understanding the excitation 
spectrum revealed in |^ is important for potential applications of the MI state. The 
zero particle number fluctuations for an ideal commensurate MI state make them 
attractive candidates for several applications, most notably as a quantum memory, 
a basis for quantum computing |20l 1211 1^ 1231 1^ I25| . and quantum simulations of 
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many-body quantum systems |25II27| . A well understood excitation spectrum can give 
valuable information about the nature and stability of the experimental approximation 
of the ideal MI state to external perturbations. 

In this paper we study the dynamics of the ID BHM under lattice modulation 
and generate excitation spectra for box and harmonic trapping of a large system over 
numerous lattice depths ranging from the SF to MI regime. We find that much of the 
features of the box system can be understood from a small exact calculation. However, 
the large system calculations were crucial for the investigation of signatures in the 
spectrum which indicate the transition from MI to SF regime in both trappings. We 
find that for a harmonically trapped system with less than unit filling the spectrum 
is similar to the commensurately filled box, but with the transition producing less 
pronounced signatures. Additional calculations progressing from the MI regime for 
a harmonically trapped system with a central filling greater than unity produces a 
spectra that has very good qualitatively agreement with the experiment. In this 
way our results are different but complementary to a very recent study of the same 
experiment by KoUath et al j28| where they discovered that excitation spectra can 
reveal information about the commensurateness of the system. 

The structure of the paper is as follows: we give an overview of the physical setup 
for the ID system in section 12.11 followed by a description of the excitation scheme 
in section 12.21 We then introduce the linear response formalism for this scheme in 
section 13.11 and in section 13.21 we give an overview of the literature describing the 
TEBD simulation method used here for the larger systems. The results are then 
presented in section^ firstly for a small box system computed exactly in section ITTl 
then for larger systems computed with the TEBD algorithm for box and harmonic 
trapping in section 14.21 and 14.31 respectively. We then end with the conclusions in 
section |S| 

2. Probing the system 

2.1. Optical lattices and the Bose-Hubbard model 

In the experiment & effective ID systems were formed from an anisotropic 3D optical 
lattice loaded with ultra-cold bosonic atoms . This is done by adiabatically exposing 
a BEC to far-off resonance standing wave laser fields in three orthogonal directions 
forming a 3D optical lattice potential VoLi'r) = J2d=i ^do sin^(te^d) where qs — 2tt/X 
and A is the wavelength of the laser light yielding a lattice period a — A/2 ,2E^. The 
height of the potential Vdo is proportional to the intensity of the d-th pair of laser 
beams, and is conveniently expressed in terms of the recoil energy Er = qg/2m for 
atoms of mass m (taking h ^ 1 throughout). 

Effective ID systems are then formed by making laser intensities in two of the 
directions r2 = y and ra = z very large (V± ~ 30Er). The confinement is then 
sufficiently strong to inhibit any tunnelling or excitations in those directions on 
the energy scales we are concerned with. The result is an array of many isolated 
effective ID systems in the ri = x direction [HI 1301 1^1^. For the remaining lattice 
intensity Vio = Vq we consider much shallower depths, but always remain deep enough 
(Vo > 4:Er) to ensure that there is an appreciable band-gap between the lowest and 
first excited Bloch band, given in the harmonic approximation by Uho = ^E^ \JVqJ E^.. 
Combined with the ultra-low temperatures of the atoms this is sufficient to ensure 
that the dynamics can be described by the lowest Bloch band of the lattice, and that 
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the tight-binding approximation is applicable [3 • 

With the centre of lattice site j in one such ID system given by Xj = ja we 
can construct a complete and orthonormal set of localized mode functions 0j(r) = 
w{x — Xj)W{y)W{z) factorized as the product of Wannier functions w and W of 
the lowest Bloch band for the shallow and deeply confined directions respectively. 
After expanding the bosonic field operator -0 (r) into these modes and restricting our 
consideration to one ID system, the resulting Hamiltonian H reduces to the ID BHM 
P] composed of M sites 

M-l M M 

where the operators hj (h-) are bosonic destruction (creation) operators for an atom 
in site j. The parameters of the BHM are functions of the lattice depth Vb with the 
matrix elements for hopping between adjacent sites j and j + 1 and on-site interaction 
strength given by ,31 

J = - j dxw{x~Xj){ -—-— + VQS\v?{qBx)\w{x-Xj+i), 



U = 2asUJho I dx \w{x ~ Xj)\ , (2) 

where a^ is the s-wave scattering length, and a Gaussian ansatz has been used for 
the tightly confined Wannier states W. The trapping offset is well approximated 
as Vj « Vrixj,!/, z), where Vrlr) describes an additional slowly varying trapping 
potential which could be due to magnetic trapping. In the case of |B] the axial 
potential of the ID system was dominated by the Gaussian beam envelopes (with 
1/e^ waists I) of the lasers for the strongly confined directions characterized by the 
trapping frequency ujt — '2Er ^ \/V±jRr PDj ■ 

The physics of the BHM is governed by the ratio U/J. Competition between 
these two terms results in a transition at temperature T = for a critical ratio 
{U/J)c from the SF to the MI regime [2131]. Mean-field theory for an infinite unit 
commensuratcly filled ID system predicts that {U/J)c ~ 2 x 5.8. However, if the 
strong quantum fluctuations present in ID are taken into account, the appearance of 
the SF regime is not predicted to occur until the critical ratio drops to {U/J)c ~ 3.85 
| 33) . The presence of trapping and the finite size of a system modifies the nature of 
the transition, prohibiting it from being sharp and so in line with the experiment [H] 
we expect the transition to occur gradually somewhere in between these limits j34| . 

2.2. Lattice modulation excitation scheme 

In the experiment j6j the ID system was initially prepared in the groundstate for 
some depth Vq, ranging from the SF to MI regime. The axial lattice depth was then 
subjected to a modulation of the form 

VoLix,t) = Vo[l + Asm{nt)]sm\qBx), (3) 

where A is the modulation amplitude as a fraction of the initial lattice depth, and 
fl is the modulation frequency |^. The modulation was applied for a fixed time 
r = 30 ms after which the energy deposited into the system was measured by time- 
of-flight imaging of central momentum width averaged over the many ID systems 
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realized. The applied modulation frequency was taken to a maximum of r2/27r = 6 kHz 
which defines the relevant energy scale for the system and was well below the band-gap. 
To compare with the experiment we take the wavelength of the light used to 
form the optical lattice as A = 826 nm, and the atomic species trapped as *^Rb, where 
Os = 5.1 nm, in all numerical values quoted. For our calculations we initially computed 
the groundstate of the system over depths U / J = 5, 6, . . . , 20 for the large systems and 
slightly shallower depths U / J = 2, 3, . . . , 20 for the small system [^, all with fixed 
particle number N . We study a small and large system with box boundaries which 
are commensurately filled as M = iV = 7 and M — N = Al respectively. We also 
consider a harmonically trapped system where Vj = muj'^x'^/2 for a slightly smaller 
system with M — 25 using ujt/^tt — 70 Hz and N = 15, as well as ojt/^tt = 100 Hz 
and N = 30. These trapping frequencies arc close to that in the experiment ED] 
where lut/2t: « 85 Hz, and where sufhcicnt to eliminate any occupation at the box 
boundaries of the system. The modulation given in equation Q was then applied 
to the system by computing time-dependent BHM parameters via equation Q. This 
includes an implicit assumption that the Wannier states describing atoms in the lattice 
adiabatically follow the variations in the lattice potential induced by the modulation. 
Given that the timescale of atomic motion in a lattice site is i^ = uJho/^Tr and that this 
is typically an order of magnitude greater than the modulation frequencies applied, the 
adiabatic assumption is reasonable. The response of the system was then measured 
via the total energy (Hq), with respects to the unperturbed BHM, for different fi. To 
demonstrate the evolution of the spectra at different depths we use the same fixed 
range for f2/27r at all depths. This range is identical to the experiment and is likewise 
quoted in kHz while the energy absorbed is expressed in units of Er- 

3. Analysis 

3.1. Linear response 

We consider a straightforward linear response treatment of this excitation scheme 
which follows under the same assumption used for the numerical calculation that 
the system is described by the BHM with time- varying parameters J[VoL{t)] and 
U[VoL{t)]- For perturbative calculations we make an additional assumption that the 
modulations are weak and approximate the variation of these functions about the 
initial depth Vq linearly resulting in a harmonic perturbation 19] 

H{t) = Ha+ AVo sm{nt) {SUHq - Jo{SJ - SU}Hj) , (4) 

where Ho = -JoHj + UoHu, SU = ^\vo, SJ = ^\vo, Uo = U[Vo] and 
Jo = J[Vo]. The perturbation is split into a part that is proportional to the 
unperturbed BHM Hamiltonian Hq and a part proportional to the hopping operator 
Hj = X^i-i (^I^j+i + h-C-) under the proviso that f/o > Jo- The first part cannot 
induce excitations and instead gives a small time-dependent shift to the unperturbed 
energies which can be ignored. As a result the excitation operator of this perturbation 
is just the hopping operator ^JE] with couphng n = AVqJq{5J — SU). By acting 
over the whole system uniformly it creates excitations with zero quasimomentum as 
expected. 

Let us label the eigenstates of the unperturbed BHM Hq as | n) with energy 
e„. The principle quantities of interest are the excitation probabilities Po^n{T,^) 
for the transitions to the excited states | n) from the groundstate | 0) due to this 
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perturbation being applied for a time r with a frequency f2. In first-order time- 
dependent perturbation theory these are given by 

F(H„(r,f]) = |4^i„|' = «:2|(n|i/,;|0)/«(r,f])|', (5) 

where In (t, il) = C dt e"^""* sin(rit) and ujna = (cn — eo)- This resuh reduces under 
the rotating wave approximation and the hmit r — ^ cxd to the famihar Golden rule 
result. We also make use of the second-order result 

2 



P^'Ur.n) 



where /i^i(T,f2) = ^^ dt J^ dt' e^'^"™*e*"'"''*' sin(m) sin(m'). The total energy 
absorbed by the system, relative to Hq, is then E{t, il) = ^„ CnPo^niT, 51) — eg. 

3.2. Numerical method 

For the large systems investigated later in sections l4 . 21 and ITTSl exact integration of the 
many-body Schrodinger equation is not feasible. To compute these results we employed 
the TEBD algorithm |3K[ 177] which is a quasi-exact numerical method that allows the 
dynamical evolution of ID quantum lattice systems with nearest-neighbor interactions 
to be computed efficiently and accurately. The algorithm has been successfully applied 
to numerous physical systems including the BHM [SHlEniEnillllllHlllS- Not long after 
being proposed by Vidal it was recognized ^9 , 43 that TEBD shares some conceptual 
and formal similarities with the well established density matrix renormalization group 
(DMRG) 2^EI method enabling the development of a new adaptive time-dependent- 
DMRG algorithm which incorporates optimizations from both. A detailed analysis of 
the accuracy and error propagation of this method was given by Gobert et al 0Tj and 
applies quite generally to TEBD also. We do not describe the TEBD algorithm here, 
except to mention some specific issues, and instead refer the reader to the relevant 
articles above for more details. 

The objective of the simulations in this paper were to map out the response of 
the system with the modulation frequency fl and depths U/J which cross the SF-MI 
transition. For sufficient sampling this required in excess of 1000 simulations for both 
the box and harmonically trapped systems presented here. Consequently we were 
limited for practical reasons to using a truncation parameter x = 30 (see |37| for 
an explanation of x) for all simulations which was lower than that strictly necessary 
in order to achieve full convergence. Thus we cannot claim that our calculations 
are quasi-exact and would not expect the many-body state given by the simulation 
at the end to have a high fidelity with the true state of the system. Instead we 
treat our calculation as an approximation in the same spirit as the Gutzwiller ansatz 
pi I4t)| (where x = 1), but with the important difference that since x > 1 we are 
permitting a non-negligible amount of quantum correlations. Given that we are only 
interested in the total energy, which is an observable composed of one- and two-particle 
correlations, we expect that this approximation should yield quantitative agreement 
for the system sizes and regimes considered here. Indeed, we have carried out more 
accurate calculations at specific points which reveal that the total energy is a robust 
observable with respects to truncation p7] . Consequently, the features of the energy 
spectrum do emerge, even with this relatively low x, and the approximations made 
do not invalidate the results presented here. For simulations of the lattice modulation 
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problem at some specific lattice depths with larger x see the recent work by Kollath 
et al EHI- 

4. Results 

4-.1. Exact calculation for a small system 

We begin by studying the exact dynamical evolution of the BHM under the lattice 
modulation for a small system composed of M = 7 sites and N = 7 atoms with box 
boundary conditions. This system is large enough to produce many of the essential 
physical features while still permitting the exact eigenstates to be computed '^. In 
fact we find that much of what is learnt from the small system can be directly applied 
to the larger systems. We make use of this by computing the energy spectrum e„ 
and perturbation matrix elements {n\Hj\0) for two lattice depths with U/J = 4 
and U/J — 20 representing the SF and MI regimes respectively. In figure QJa) the 
MI spectrum is shown with its characteristic gapped structure composed of Hubbard 
bands located around multiplies of the dominant interaction energy U and spread 
by finite hopping J. In addition to the spectrum, the perturbation matrix elements 
{n\Hj\0) which are of order 0{J/U) or above are shown as the vertical lines. For 
the MI the most numerous (and strongest) contributions are to the [/-Hubbard 
band which is described by 1-particle-hole (1-ph) excitations like those depicted in 
figure ^c)(i). The matrix elements to the 2J7-Hubbard band, which is composed 
of two 1-ph excitations shown in figure ^c)(ii), cancel to first order explaining the 
absence of lines for this manifold. However, there are a small number of matrix 
elements to first order connecting the groundstate to the 3C/-Hubbard band via two 
1-ph excitations with both particles on the same site as in figure ^c)(iii), but not to 
three 1-ph excitations as in figure nic)(iv). In figure ^b) the 'gapless' SF spectrum 
is shown. Here in contrast to the MI there are significant contributions to {n \Hj\ 0} 
stretching from below an energy of U to below 3C/. Separated from this there are 
contributions tightly distributed around 4C/. Given the comparatively equal strengths 
of the hopping and interaction terms for the SF regime at [// J = 4 there is no simple 
picture of either the groundstate or the excitations related to these contributions as 
there was for the MI regime. Despite this the exact eigenstates for this small system do 
reveal two important details: firstly two 1-ph configurations like those in figurenic)(iii) 
have an average energy relative to the U/J = 4 groundstate which exceeds AU; and 
secondly these types of configurations are the dominant contributions in the relevant 
eigenstates around 4C/. This overlap in the nature of the excitations points to the 
possibility, which we shall shortly confirm, that a resonance to 3C/ excitations in the 
MI regime will evolve into a 4C/ resonance as the SF regime is entered. 

We first consider the modulation scheme with a weak amplitude A = 0.01 where 
linear-response is applicable. In figure 12 (a) the total energy of the system e is shown 
after the modulation as a function of fl for the MI regime with U/J = 20. As 
expected from the matrix elements (n|_ff,/|0} in figure ^^a) there is strong response 
centred around U spread by the width of the t/- Hubbard band which is approximately 
10 J ~ 0.8 [kHz], and a smaller (barely visible) response at 3U. The linear response 
predictions based on equation Q are shown also and agree well. For the largest 
peak the growth of energy over the modulation time is shown in figure |2Ib) and the 
slight overestimation of the energy by linear-response is evident at later times. The 
same calculation for the SF regime with U/J = 4 is shown in figure |2Ib). Again 
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Figure 1. The energy spectra for the first 300 eigenstates n for (a) U/J = 20 MI 
regime, and (b) U/J = 4 SF regime. In both cases the vertical (red) lines denote 
the presence of a matrix element | (n \Hj\ 0) | connecting an excited state to the 
groundstate which is of order 0{J/U) or above. In (c) a schematic depiction of 
the excitations which exist in the (i) U-, (ii) 2U- and (iii-iv) 3C/-Hubbard bands 
of the MI regime in (a) is given. 



the results of linear-response agree well and the structure follows that of perturbation 
matrix elements with a broader response between U and SU and an equivalent strength 
resonance at AU. 

In line with the experiment |S] from this point on we consider much stronger 
modulations with A = 0.2. The result of these for the MI and SF regimes is given 
in figure |3| For the MI in figure Ufa) we see that the discrete resonances around 
U have now filled out into a single peak centred on U that is about 25% wider 
than the [/-Hubbard band at ~ 1 [kHz], and that the stronger modulation has now 
increased height of the 3U peak relative to the U peak. In figure El^a) we have also 
included both first- and second-order perturbation theory results. For such strong 
modulations the applicability of perturbation theory is highly questionable, especially 
for the long times considered here. This is exemplified by the gross overestimation of 
the central peak at U by both hnear- and quadratic-response. In figure Ol^b) curve 
(i) shows the saturation of energy absorption for the central peak at U and departure 
from linear response after a short time. In contrast for the peak at 3U figure |3Ib) 
curve (ii) shows that linear-response underestimates the energy absorbed due to its 
neglect of the role these eigenstates play in the indirect processes to higher energies. 
However, the use of quadratic-response given by equation ^ provides some useful 
insight into the additional structure seen for the strong modulation which do not 
appear for weak modulations and is not predicted by linear-response. Specifically, the 
small satellite peaks either side of the dominant peak at U, which are located at U/2 
and 3U/2 respectively, only appear at second and higher order leading to their natural 
identification as the absorption of two quanta of energy Q from the perturbation to 
reach to U- and 3t/-Hubbard bands. We note also that even under strong modulations 
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Figure 2. For a weak modulation of the small system we have (a) the total energy 
e absorbed in the MI regime with U / J = 20 as a function of f2, (b) e over time for 
the largest MI peak at f2/27r ~ 1.7 [kHz] for the exact (E) and linear-response (L) 
calculations, and (c) e against O for the SF regime with U / J = 4. In (a) and (c) 
the markers (o) and (*) denote the exact and linear response results respectively, 
while the lines are drawn to guide the eye. 



no resonance at 2t/ is seen in agreement with recent findings j28| for commensurate 
systems. For the SF regime in figure El^c) a broad response is seen spanning the 
region between U and 3C/ consistent with linear response, but higher-order effects 
have resulted in saturation and merging of linear response peaks. As with weaker 
modulations there is a separate resonance centred at 4C/ which is now as equivalent 
in strength and broader. 

With the full exact calculation we have computed the response of the system for a 
sequence of depths C// J = 2, 3, . . . , 20 which are displayed in figure01(a) demonstrating 
the evolution of the MI spectrum in figure |31(a) into SF spectrum in figure El (c) with 
decreasing lattice depth. By examining the accompanying colour-map in figure 0](b) 
we can make a number of initial observations regarding the changing characteristics 
of the spectrum over the SF-MI transition. Firstly, as with the weak modulations, the 
MI peak at U is seen to broaden and shift upwards in energy into the SF response, and 
this change is most prominent after the depth is lower than U / J ^ 12. Secondly, it is 
now clear that the resonance at 3[/ for deep in the MI regime reduces only slightly in 
energy and ends up as the 4U resonance in the SF regime. Additionally, in line with 
weaker modulations in figure El the 4U peak is stronger than the 3C/ peak it evolves 
from, and signatures of this change are already visible in figure 01 (a) . Finally, by 
progressing to shghtly shallower depths {U/ J = 2, 3) in figure^we see the broadening 
of the SF spectrum over the entire frequency range considered caused by the eventual 
merging of the U to 3C/ response with the 4C/ peak. 
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Figure 3. For a strong modulation of the small system we have (a) the total 
energy e absorbed in the MI regime with V j J = 20 as a function of O, (b) e over 
the time at the MI peaks (i) V and (ii) 3(7 for the exact (E) and linear-response 
(L) calculations, and (c) e against f2 in the SF regime with U j J = 4. For (a) and 
(c) the markers (o), {*) denote the exact and linear-response results respectively, 
while in (a) (O) denotes the quadratic-response. The lines arc drawn to guide the 
eye. 
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Figure 4. (a) The total energy e absorbed by a small box system over a sequence 
of depths, ranging from the MI to SF regimes, as a function of SI for strong 
modulations, (b) A colour-map of plot (a). 



4-2. Large system in a box 

To put some of the observations made in the previous section on a firmer footing we 
now examine a larger system with M = 41 and N ~ 41, again with box boundaries, 
under the same modulation scheme with A = 0.2. To compute the dynamical evolution 
for this system we resort to the TEBD algorithm using a lattice site dimension rf^ = 5. 
In figure |31 the one-particle density matrix pjk = {b^jbk) for the groundstate of (a) 
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the MI with U/J = 20 and (b) the SF with U/J = 5 is plotted, along with the 
corresponding spectra obtained from these groundstates in (c) and (d) respectively. 
As might be expected the increase in the size of the system reduces finite size effects 
resulting in much smoother excitation profiles. However, the essential features of 
these plots still follow directly from our analysis of a small system in terms of both 
the position and width of the resonances. Specifically, the width of the U peak in 
figure inic) is again slightly larger than the width of the [/-Hubbard band for this 
system given approximately by first order perturbation theory as 12 J ^ 1 [kHz]. The 
effects of saturation also appear to follow similarly with the maximum energy absorbed 
per particle being nearly identical for the two system sizes. 



\pjk\ fS 




Figure 5. The one-particle density matrix pjfc of the groundstate for (a) the MI 
regime with U/J = 20 and (b) the SF regime with U/J = 5. The corresponding 
excitation spectra of the total energy e against C for these groundstate is shown 
in (c) and (d) respectively. 



The absorption spectrum over a sequence of depths U/J = 5, 6, . . . , 20 ranging 
from the SF to MI regime are shown in figure El The increased size of the system 
permits us to investigate more conclusively some of the observations made earlier 
regarding the evolution of the absorption spectra with decreasing depth. Firstly, we 
focus on the maximum strength of the 3C/ -^ 4X1 peak with U/J which is shown in 
figure[3(a). The plot shows that the strength of this peak displays a distinct alteration 
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in its behavior close to U / J ~ 10 where it jumps from an increasing to a decreasing 
curve. This behaviour is consistent with what was seen earher for the small system 
and can be attributed to the increasing SF character with a strong AU peak becoming 
dominant over the weaker 3C/ MI peak. 
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Figure 6. (a) The total energy e absorbed by a large box system over a sequence 
of depths, ranging from the MI to SF regimes, as a function of fl for strong 
modulations, (b) A colour- map of plot (a). 
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Figure 7. (a) The maximum energy e absorbed for the 3{/ — > 4f7 peak with 
the lattice depth U / J . (b) The centre position He, and (c) width Om of the U 
resonance with the lattice depth U /J . The dotted lines in all cases is drawn to 
guide the eye. 



To understand the changes exhibited by the resonance around U we fit this peak 
with a smoothed-box function of the form 



e{n) 



-M- 



cmax 



(7) 



l_)_e[(n-n,)2-«,2]/s2 ' 

where Vic, w and s are parameters specifying the center, top width and step size of 
the profile respectively, whilst A^ = 1 + e^™ Z'* is the scaling factor. The advantages 
of this function (which are even more apparent in the next section) is that it faithfully 
describes both the width and centre location of broad topped resonances. The centre 
fie and the full- width-half-maximum fi^, extracted from this fitting then characterizes 
quantitatively the behaviour which is evident in figure El(b). 
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In figure[7|(b) we see that the centre Vt^ of the U peak remains relatively stationary 
until U / J ^ 12 where it then shifts more dramatically to higher energies. This is again 
a signature of increasing SF character consistent with the disappearance of the gap 
between the U- and 2C/-Hubbard bands and the shifting of the strongest contributions 
to (n \Hj\ 0) to higher frequency seen for the small system. The width fi^ in figure[3(c) 
instead displays a gradual increase without any pronounced changes. 

It is clear from the structure of both the MI and SF spectra shown in figure Efc) 
and (d) that the box system has some differences compared to the ID spectra obtained 
in the experiment. These differences are not surprising given that the SF state of 
the box system at U / J — 5 exhibits significant quantum depletion of 75%, which 
is much larger than that in the experiment. Also, at all depths considered here the 
groundstate of the box has a near homogeneous density of 1 atom per site, as can be 
seen in figure El^a) and (b), which is different from the harmonically trapped system 
in the experiment. 

^.3. Large harmonically trapped system 

To approach a setup closer to that of the experiment we now consider a system with 
M — 25 and A'^ = 15 superimposed with harmonic trapping using loj^ J^-k = 70 Hz. A 
significant difference is that at U / J = 5 the depletion of the SF state is reduced to 
50% in line with the experiment for this setup. The one-particle density matrices for 
the Ml and SF groundstates in this case are shown in figure|SJa) and (b) respectively. 
We now see that the MI state is composed of a central core with unit filling along 
with small SF lobes at the edges. This change in structure compared to figure Efa) 
introduces additional types of excitations such as those caused by particles hopping 
into the vacuum surrounding the MI core. For the SF state there are now significantly 
greater off-diagonal correlations in pjk compared to the box system. Despite these 
differences, however, the form of the spectrum shown in figure ISfc) and (d) remains 
very close to the box system displaying the same characteristic peaks as before. The 
effects of trapping manifest themselves in this case by flattening of the U peak, along 
with the broadening and shifting of the AU peak. The maximum energy per particle 
again remains approximately equal to the previous cases considered, showing that the 
saturation effects are not significantly altered by the trapping. 

In figure El we report the spectra for the system ranging over SF to Ml regimes 
as before. Again, the essential features of this plot follow from our earlier discussion. 
Performing the same analysis on the evolution of the spectrum demonstrates that 
maximum strength of the 3C/ -^ AU peak has a maximum at U / J ^^ 8 as shown in 
figure [TUT a). The trapping is also seen to introduce peaks and troughs to this curve 
compared to the smooth behaviour of the box system. Also, in contrast to the box 
both the centre f2c and width ^^ of the U resonance show a gradual change from the 
MI to SF character in figure lTnT b) and (c). This refiects a difference in the evolution of 
U resonance, which for the box develops a visible two-peaked structure in figurej^a), 
whereas in the harmonic trap this is washed out into a broader and fiatter response 
in figure 13 a). 

We have seen that the introduction of harmonic trapping for the case considered 
above has not lead to any new features in the spectrum beyond a small amount of 
shifting and broadening. One such feature which might be expected is the emergence 
of a resonance at 2f/. However, it is known that the relevant excitation processes for 
a 2U resonance require either excitations to already be present in the system due to 
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Figure 8. The one-particle density matrix pji^ of the groundstate for a trapped 
system with o^t/Stt = 70 Hz for (a) the MI regime with U/J = 20 and (b) the SF 
regime with U/J = 5. The corresponding excitation spectra of the total energy 
e against the modulation frequency Q for these groundstate are shown in (c) and 
(d) respectively. 



finite temperature |13| . or significant inhomogeneity in the density caused by a trap 
|28|. Since the system above is at T = and has virtually no incommensurability in the 
MI, beyond the small SF lobes at the edges, the lack of a 2U resonance is consistent. 
In the experiment both finite-temperature and inhomogeneity are expected to make 
important contributions to the spectrum. Presently there are still open questions as to 
which is the most dominant and what interplay there might be between these effects. 
Future studies with the TEBD algorithm, exploiting recent developments which permit 
the simulation of master equation evolution of ID finite-temperature and dissipative 
systems 0^1 j could allow both these effects to be treated on an equal footing. 

To demonstrate the contribution of inhomogeneity we have performed additional 
calculations again using Af — 25 but with a larger occupancy of A^ = 30 and larger 
trapping frequency ujt/2tt — 100 Hz. In ffgure irTT a) the one-particle density matrix of 
the MI groundstate at U/J = 20 is shown and displays the coexistence of significant 
MI and SF regions characteristic of trapped systems. Despite using a larger lattice site 
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Figure 9. (a) The total energy t absorbed by a M = 25 and N = lb harmonically 
trapped system over a sequence of depths, ranging from the MI to SF regimes, as 
a function of C for strong modulations, (b) A colour-map of plot (a). 




Figure 10. (a) The maximum energy e absorbed for the 31/ — > 4{/ peak with 
the lattice depth U / J . (b) The centre position Vic, and (c) width O^ of the U 
resonance with the lattice depth U /J . The dotted lines in all cases is drawn to 
guide the eye. 



dimension dg — Q for this setup we found that truncation in on-site occupancy Hmited 
the smallest U / J from which we could reliably compute the spectrum to U / J — 9. 
The one-particle density matrix for the U I J = 9 groundstate is shown in figure [TTT b') 
and appears to be dominated by a central SF region. The excitation spectra for 
these two groundstates are shown in figure Hir e') and (d) respectively. For the MI 
in figure lllf c) a pronounced 2U resonance is seen in addition to the peaks seen in 
previous MI spectra. The appearance of a 2U peak in the T — Q case examined here 
arises predominantly from particles in the unit filled MI shell hopping into the doubly 
occupied region at the centre 28 , 6 which is evident in the one-particle density matrix 
shown in figure ITIT a'). For the shallower depth U/J = 9 the spectrum in figure [TlT d) 
exhibits the initial merging of the U and AU peaks with this 2U peak. In figure IT^ the 
evolution of the spectrum between these limits is shown. The form of these spectra 
is very reminiscent of that shown in the experiment [HI and it is likely that a more 
detailed analysis of this setup to shallower lattice depths would reveal the broad SF 
resonance seen experimentally. 
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Figure 11. The one-particle density matrix pj^ of the groundstate for a trapped 
system with ojt/Stt = 100 Hz for (a) the MI regime with U/J = 20 and (b) 
an intermediate regime with U/J = 9. The corresponding excitation spectra of 
the total energy e against the modulation frequency H for these groundstate are 
shown in (c) and (d) respectively. 



5. Conclusions 



We have examined in detail the dynamical response at T = of the ultracold atoms in 
an optical lattice subjected to lattice modulations. We have reported the evolution of 
the excitation spectrum with the lattice depth from the SF to MI regime for small and 
large box systems, and for large harmonically trapped systems. For the box system 
we identified two pronounced signatures of the transition from the MI to SF regime, 
specifically the strength of the iU —^ 411 peak and the centre of the U resonance. The 
introduction of trapping for the case considered here, where the occupancy remains 
less than or equal to unity, does not significantly alter the evolution of the spectrum 
compared to a box system. While we find that it does wash out any pronounced 
changes in the structure of the U resonance with depth, the strength of the 3U -^ 411 
peak still exhibits a signature of increasing SF character. 

We also presented calculations showing spectra progressing from the MI regime 
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Figure 12. (a) The total energy e absorbed by a Af = 25, A' = 30 harmonically 
trapped system over a sequence of depths, ranging from the MI to SF regimes, as 
a function of Q for strong modulations, (b) A colour-map of plot (a). 



for a harmonically trapped system which has a central region with greater than unit 
filling. We found that this changes the structure of the spectrum by introducing a 
2U resonance and brings our results closer with the experiment. Direct comparison 
to the experiment fS-, however, is difhcult due to the 3D rethermalization performed 
prior to measurement, which cannot be simulated with TEBD, and the fact that the 
measurement itself was an averaged result over many ID systems in parallel with 
differing total particle number. Despite this our results point to that fact that the 
BHM is sufficient to explain all the features discovered in the experiment and that the 
experiment was a clean realization of this model as expected. 
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